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Abstract 

Parametric instabilities of parallel propagating, circularly polarized finite amplitude Alfven waves 
in a uniform background plasma is studied, within a framework of one-dimensional Vlasov descrip- 
tion for ions and massless electron fluid, so that kinetic perturbations in the longitudinal direction 
(ion Landau damping) are included. The present formulation also includes the Hall effect. The 
obtained results agree well with relevant analysis in the past, suggesting that kinetic effects in the 
longitudinal direction play essential roles in the parametric instabilities of Alfven waves when the 
kinetic effects react "passively". Furthermore, existence of the kinetic parametric instabilities is 
confirmed for the regime with small wave number daughter waves. Growth rates of these insta- 
bilities are sensitive to ion temperature. The formulation and results demonstrated here can be 
applied to Alfven waves observed in the solar wind and in the earth's foreshock region, (accepted 
to Physics of Plasmas) 
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Circularly polarized, finite amplitude Alfven waves are ubiquitous in collisionless space 
plasmas, for example in the solar wind, and also in the foreshock region of planetary 
bowshocks where ions backstreaming from the bowshock 1 generate large amplitude Alfven 
waves 2 ' 3 . Parametric instabilities are of particular interest for dissipation of quasi-parallel 
Alfven waves there, since they are typically robust for linear ion-cyclotron damping (due 
to small wave frequencies) and for linear Landau damping (due to small propagation angle 
relative to the background magnetic field). 

Within a framework of the Hall-magnetohydrodynamic (MHD) equations, a number of 
research has been carried out on the parametric instabilities of parallel propagating, circu- 
larly polarized finite amplitude Alfven waves 4-7 , which is an exact solution of the nonlinear 
Hall-MHD equation sets. Among various parametric instabilities, modulational, decay, and 
the beat instabilities are important. In general, growth rates of the instabilities depend on 
the parent wave amplitude, wave number, polarization, and the ratio of the sound to the 
Alfven velocity (= C S /C A = y/fr). 

While a large number of past studies deal with the instabilities within the framework of 
fluid theory, since the ion beta (the ion fluid to magnetic pressure ratio) is not small in the 
solar wind and in the earth's foreshock region, the ion kinetic effects cannot be neglected. 
From this perspective, numerical analysis using the hybrid simulation code (super-particle 
ions + massless electron fluid) was performed by Terasawa et. al. 7 and Vasquez 8 . Mjolhus 
and Wyller 9 ' 10 derived and discussed in detail a nonlinear evolution equation (kinetically 
modified derivative nonlinear Schrodinger equation, first derived by Rogister 11 ), which de- 
scribes the wave modulation of weakly despersive, nonlinear Alfven waves including the ion 
Landau damping. Fla et. al. 12 and Spangler 13 ' 14 discussed the modulational instability of 
parallel propagating, circularly polarized Alfven waves, and pointed out that, while the ki- 
netic effects suppressed the parametric instabilities present in the fluid model ( "conservative 
modulational instability (CMI)"), the ion kinetics could evoke a new instability ("resonant 
particle modulational instability (RPMI)") 12 . Inhester 15 first showed the kinetic formulation 
for parametric decay instabilities without dispersion. Gomberoff 16 and Araneda 17 discussed 
the ion kinetic effects on parametric instabilities phenomenologically. Common view among 
past analytical works on parametric instabilities is that the ion kinetic effects expand the 
unstable parameter regime and reduce the maximum growth rate. Recently, Passot and 
Sulem 18 derived the dispersive Landau fluid model, and made comparison with some of the 
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other models in the past from both analytical and numerical standpoints 18 ' 19 . 

Linear analysis play important role in verifying the instabilities seen in the numerical 
simulation 8 and provide a more complete reporting of growth rates as a function of the 
parameters. To our knowledge, linear kinetic analysis with dispersion has not yet been 
demonstrated. Our aim in this study is to derive a kinetic formulation including dispersion 
systematically, and present a linear analysis and compare it with other analytical results. 
Our treatment of the ion kinetic effects is similar to that in Munoz et. a/. 20 , in which 
relativistic dispersion relation for parametric instabilities in an electron-positron plasma is 
derived. 

We consider the parametric instabilities of parallel propagating, circularly polarized finite 
amplitude Alfven waves, which are the exact solution within the Hall-MHD equation sets. 
Assuming weak ion cyclotron damping, we include the kinetic effects only along the longitu- 
dinal (x) direction. Let f(x,t,v) be the ion distribution function, and define the integrated 
longitudinal distribution function, 

g(x,t,v x ) = J f(x,t,v)dv y dv z . (1) 

Then the governing set of equations can be written as, 

du du 1 db 

dt x dx pdx ) 

db d I ~ _ % db\ 

at = -dx-l Uxb - u + -pdx-)' (3) 
m = -^-( e * + V^A)^, (4) 

1 d f\b\ 2 \ 

e x = -u y b z + u z b y - -— ( — + p e \ , (5) 

where p is the plasma density (quasi-neutrality assumed), u = (u x ,u y ,u z ) is the ion bulk 
velocity vector, b = b y + ib z and u = u y + iu z are the complex transverse magnetic field and 
bulk velocity, respectively, and e x is the longitudinal electric field. All the normalizations 
have been made using the background constant magnetic field, density, Alfven velocity, and 
the ion gyro-frequency defined at a certain reference point. 

The total pressure is given as p = p e + Pi, where isothermal electrons are assumed, i.e., 
p e = T e p. Also it is assumed that the ion and electron pressures are isotropic both at 
the zeroth and at the perturbation orders. Usual beta ratio slightly differs from the ratio 



between the sound and the Alfven wave speed squared, = (jePe + 7iA)/2, where j e — 1 
and 7i are the ratios of specific heats for electrons and ions, respectively. 

At the zeroth order we consider the parallel propagating Alfven wave given as 

b p = b exp(i(u) t- k x)), (6) 
u p = u exp(i(uj t - k Q x)) , (7) 

with u = —bo/v^o, phase velocity v^o = ooo/k , p = 1, u x o = 0, together with the zeroth 
order dispersion relation 

ul = kld+uv). (8) 

We adopt the notation that the positive (negative) lu corresponds to the right- (left-) hand 
polarized waves. For the zeroth order longitudinal distribution function, we assume go(v x ) = 
expi-vl/vU/y/nvth {y 2 th = 

Then we add small fluctuations given by 

5E = E + exp(i(f) + ) + E_ exp(i0_), (9) 
SF = -(Fiexp(«$) +c.c), (10) 

where E represents the transverse variables (b and u), and F represents the longitudinal 
variables (g, p, and u x ), <f>± = uj±t — k±x, $ = Qt — Kx, k± — ko ± K, uj± = uo ± fi, and 
c.c. denotes the complex conjugate. 

Then, (2)- (5) at the first order produce 

(u)± -k 2 ±- u±k±)b^ = ^(b u± - u k )u xl 

b k k± 

— (w± + l)pi, (11) 

K[b (b + + b*_) + T ePl }dg (v x ) 
i£ — v x K ov x 

where T e = (3 e /2 and the asterisk denotes the complex conjugate. 
Integration of the above yields 

Pl = -K[b (b + + b*_)+T ePl ]D, 
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where £ = Q/v t h/K, and Z(£) is the plasma dispersion function. This equation can be 
written as 

Pl= A (K+i >_), A = - T ^ 5 . (14) 

In a similar way, we have 

u xl = B(b+ + b*_), B = -K(b + T e A)C, (15) 

where n = 1. Combing (11)-(16), we obtain 

= L + P_ + L_P+ (17) 

where 





= ul-k%(l+u±), 


(18) 


s± 


b k± B b k k± 
2 2 


(19) 


T± 


Vok k± b k k± 

= d H A, 

2 2 


(20) 


P± 


= S±L±-T±. 


(21) 



By taking the cold limit ~ 0, i.e., £ >> 1), the dispersion relation of the Alfven wave 
parametric instabilities in the Hall-MHD is obtained 7 . We note that (17) is an odd function 
with respect to both Q and K. 

Now we examine the numerical solutions of (17). First, we refer to the results by 
Inhester 15 , in which the kinetic dispersion relation of parametric decay instability is de- 
rived without dispersion, using drift kinetic model. In order to make comparison with their 
results, we omit the dispersion terms in (17): the third term in the RHS of (18) and the 
second term in the RHS of (19). Note that in the dispersionless system, the modulational 
instability (both CMI and RPMI) does not occur and the growth rate of the decay instability 
differ from one in the dispersive system 7 . By further assuming the cold limit, our equation 
leads to the MHD dispersion relation obtained by Goldstein 21 and Derby 22 , (but not to the 
model by Inhester (1990)). On the other hand, numerically obtained dispersion relations 
based on our model and that of Inhester at least qualitatively agrees, in the sense that the 
ion kinetic effects enlarges the unstable parameter regime and also the maximum growth 
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FIG. 1: Growth rates of the density modes (whose frequency /wave number is normalized to the 
parent wave parameters), driven by the decay instability of dispersionless Alfven waves with am- 
plitude b = 0.447, = 0.6, and T r (= Ti/T e ) = (a) (=fluid model, black cross), (b) 0.03 (blue 
broken line), (c) 0.2 (light blue), (d) 1 (green dotted line), (e) 5 (orange), (f) oo(T e = 0) (red solid 
line). The results shown here are in agreement with Inhester (1990) and Bugnon et. al. (2004). 

rates are reduced (Figure 1). These results are also in qualitative agreement with some of 
the past works 15 ' 19 . 

We now turn our attention to the parametric instabilities of the dispersive Alfven waves. 
Figure 2 shows the growth rates computed from (17) for right-hand polarized Alfven waves, 
for various temperature ratios, T r (— Ti/T e ). The decay-like instabilities have positive growth 
rates when K > \Jcq\ = 0.408. The results here qualitatively agrees with numerical results 
obtained by Vasquez 8 and Bugnon 19 (Table 1 in both papers), which suggests that the ion 
kinetic effects in the longitudinal direction play essential roles in kinetic modification of the 
Alfven parametric instabilities, when $ is relatively small, i.e., the kinetic effects "passively" 
react to the fluid dynamics. Figure 3 shows the growth rates under the same parameters as 
in Figure 2 except that k = 0.102 and u = 0.107 (R-mode). Both regimes K > \k \ and 
K < |fco| are plotted. The former corresponds to the "conservative decay instability (GDI)", 
which has the finite growth rate at T r = as we see in Figures 1 and 2. On the other hand, 
the latter, the RPMI, is destabilized only for finite T r , although the growth rate is typically 
1 or 2 orders less than the CMI. The RPMI is quenched when T r = 0, suggesting that the 
instability is a product of Landau resonant effects (Fla et. al. 12 , Spangler 13 ' 14 ). The RPMI 
exists even at small wavenumbers. 

As shown in Figure 3 (c)-(g), the growth rate of the RPMI is sensitive to /3j. In particular, 
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FIG. 2: Growth rates of the density modes (frequency/wave number is normalized to the ion- 
gyrofrequency/ion inertial length), driven by the decay instability of dispersive Alfven waves with 
6 = 0.5, k = 0.408, ojq = 0.5 (R-mode), = 0.45, and T r {= Ti/T e ) = (a) (=fluid model, 
black cross), (b) 0.0227 (blue broken line), (c) 0.36 (light blue), (d) 1 (green dotted line), (e) 2.75 
(orange), (f) oo(T e = 0)(red solid line). The results agree well with Vasquez (1995) and Bugnon 
(2004). 
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FIG. 3: Same as Figure 2 except that b = 0.5, k = 0.102, u = 0.107 (R-mode), j3 = 0.45, with 
T r (= Ti/T e ) = (a) (=fluid model, black cross), (b) 0.0227 (blue broken line), (c) l(f3 e = Pi = 
0.225) (green dotted line), (d) oo(T e = 0) (red solid line), and when fa = 0.225 with (e) fa = 0.5 
(square), (f) fa = 1.0 (circle), (g)fa = 2.0 (diamond), respectively. 

as Pi is increased, the growth rate of the RPMI is enhanced while that of the CDI is reduced. 
When Pi = 2 and p e = 0.225, the growth rates of both instabilities become comparable. The 
CDI growth rate is increased as ftf (=/3 e /2 for the present case) is reduced. The decrease of 
the growth rates of the CD(M)I is relaxed at large T r (c.f., Fig.l, 2 (e), (f)). 

Finally, we discuss briefly the instability of left-hand polarized Alfven waves (Figure 
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FIG. 4: Same as Figure 2 except that bo = 0.4, ko = 0.56, u>o = —0.425 (L-mode), (3 = 0.6, and 
T r (= Ti/T e ) = (a) (=fluid model, black cross), (b) 0.03 (blue broken line), (c) 0.2 (light blue), 
(d) 1 (green dotted line), (e) 5 (orange), and (f) oo(T e = 0) (red solid line). 

4). Under the parameters used, the CMI (the maximum growth rate ~ 0.4), the CDI 
(likewise, ~ 0.6), and the beat instability (likewise, ~ 0.8) are driven unstable at T r = 0. 
Figure 4 shows that the growth rate of the decay-like instability (~ 0.8) exceeds that of the 
modulational-like instability at T r >~ 0.2. 

In the present paper we discussed the kinetically modified parametric instabilities of cir- 
cularly polarized, parallel-propagating Alfven waves within a framework of one-dimensional 
system, which includes longitudinal kinetic perturbations. The obtained dispersion relation 
(17) is numerically evaluated, and compared with relevant works in the literature. The real 
part of Q in (17) is not shown here, but we have confirmed that the propagation direction of 
the sideband Alfven waves and the density fluctuations are consistent with the past studies. 

While our analysis includes $ as a free parameter, it should be retained at a modest 
value, since at high kinetic response of the transverse distribution function cannot be 
neglected, and the kinetic effects must be included in the parent wave dispersion relation 
(9) 8 ' 19 ' 23 ' 24 . Comparison of Fig.2 with past studies suggests that our assumption on basic 
equations is valid when = 0.45, but is not for j3i ~> 1.5 8 . This remark also applies 
to simulation studies using the hybrid code: it is important to use kinetic instead of fluid 
dispersion relation to give initial wave when the ion beta is large. 

Finally we point out that, despite the importance of the parametric instabilities in the so- 
lar wind or in the earth's foreshock, their presence has not yet been clearly demonstrated by 
spacecraft experiments. Linear theory predicts that parallel propagating waves are mainly 
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excited in the solar wind or in the earth's foreshock 25 ' 26 . However, recent multi-point mea- 
surement observations suggest that these waves tend to propagate obliquely to the back- 
ground magnetic field 1 ' 27 ' 28 . Numerical simulation study using hybrid code 3 should be effec- 
tive in understanding the parametric instabilities of obliquely propagating large amplitude 
MHD waves. We hope the situation to change via knowledge of kinetic modification of the 
instabilities as well as state-of-the-art data acquisition and analysis techniques. 

We thank Drs. S. Matsukiyo, V. Munoz, and T. Passot for fruitful discussions and com- 
ments. This paper has been supported by JSPS Research Fellowships for Young Scientists 
in Japan. 
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